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SUMMARY 

After a review of existing publications on wind-tunnel correct- 
ions at high Mach numbers, an approximate method is given for deter- 
mining the corrections due to model displacement and wake displacement 
behind resistance bodies and due to the lift. The correction compu- 
tations are first carried out for the incompressible flow. According 
to the Prandtl principle, the models and the wind tunnels in the 
incompressible flow are made to correspond to the models and the wind 
tunnels in the compressible flow,* between the corrections of these two 
correlated tunnels definite relations then exist. Because of the Prandtl 
principle is applied only to the flow at a large distance from the 
models, the wind-tunnel corrections can also be computed if the 
assumptions in the Prandtl principle are not satisfied in the 
neighborhood of the model . The relations are investigated with 
particular detail for fuselages and wings of various spans in closed 
circular tunnels . At the end of the report a comparison is made 
between the computations and the tests in the DVL high-speed wind 
tunnel . 


I . STATEMENT OF PROBLEM AND REVIEW OF LITERATURE 

If a model of finite thickness Is mounted in a closed tunnel, 
the tunnel at the model cross section is narrowed by a definite amount 
because of the displacement of the model. The air is now forced by 
this contraction to flow around the model with greater velocity than 
would be the case in an unlimited air stream. . In the application of 
the results of measurement to the free air stream, the measured values 
are therefore to be correlated with a higher velocity than corresponds 
to the velocity in the tunnel without the model . 

+ "Windkanalkorrekturen bei hohen Unterschallgeschwindigkeiten unter 
besonderer Berucksichtigung des geschlossenen Kreiskan^ls . " 
Forschungsbericht Nr. 1216, Deutsche Versuchsanstalt fur Luftfahrt, 

E. V., Institut fflr Aerodynamik, Berlin-Adlershof , May, 1940. 
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For wind tunnels with small velocities, as have so far been 
predominantly used, the increase of the velocity due to the model is 
so small in view of the generally small model dimensions that it may 
generally be neglected. For higher tunnel velocities the velocity 
correction increases very rapidly, however, and at a more rapid rate 
the more nearly the velocity approaches that of sound. This fun dam ent al 
behavior is shown, for example, by Ferri (reference 1, p. 112) by 
computing the velocity increase due to a definite amount of tunnel 
narrowing on the assumption of a uniform velocity distribution over 
the narrowed cross section (see fig. l). From these computations it 
is found, for example, that a narrowing by the model of the tunnel by 
1 percent gives, at small Mach numbers, an increase in velocity of 
1 percent^ for the Mach number M = 0.8, 2.7 percent,- and for M = 0.90, 
as much as 11.0 percent. 

This simple rough computation indicates the absolute necessity, 
on the one hand, of computations for the correction of the tunnel 
velocity due to the model displacement and, on the other hand, the need 
for choosing the dimensions of the model in relation to those of the 
tunnel so as to be smaller the nearer the velocity of sound is 
approached . 

Accurate computations for flows in a compressible fluid, on 
account of the complicated considerations, currently offer little 
promise of success when it is considered that even the simple case of 
the two-dimensional flow about a wing in an infinite compressible air 
stream already lies at the limit of present-day computation possibili- 
ties . Approximations must therefore be sought that come as close ' 
as possible to the actual facts. 


1 . Approximate Computation of Lamia 

A simple estimate of the correction for two-dimensional flow was 
given by Lamia (reference 2), who in his computations took into 
account the compressibility of the air up to terms of higher order and 
made use of the following concept as a basis : 

If a wing with infinite span is placed in an air stream (without 
lift), two streamlines lying symmetrical to the wing at a large distance 
ahead of and behind the wing are separated by a distance bop, while in 
the plane of the wing they are displaced by the greater distance 
hoo + Ah. Between two streamlines at distance h^ therefore, the 
flow in the model cross section no longer transports the same quantity 
of air Go, as between the streamlines the same distance apart at a 
large distance ahead of or behind the wing, but a smaller quantity of 
air Gob - AG. Lamia takes the approach velocity v in the free air 
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stream to be increased by the amount Avq until at the model cross 
section between two streamlines at distance h^ the quantity G*, 
is transported. The required increase in the velocity is considered 
as an approximate value for the velocity correction to be applied to 
the closed tunnel; the distance between the two streamlines h® 
are set equal to the distance between the tunnel walls. 

An important result of this estimate of Lamia was that the 
velocity corrections given by Ferri were recognized as upper limiting 
values, which can be attained only for very slender models with very 
large chord in comparison with the tunnel diameter. In by far the 
majority of cases, however, the chord of the model is so small in 
comparison with the tunnel diameter and its thickness ratio so large 
that the disturbance velocities produced by the model from the 
immediate neighborhood of the model up to the tunnel wall are greatly 
decreased and therefore also the effect on the flow produced by the 
wall becomes considerably smaller. For example, with 4 percent cross- 
sectional narrowing by the model, the velocity correction for the Mach 
number M = 0.75 is given by (from fig. 8 of reference 2) : 


Model 

Thickness ratio, 
d/t 

Velocity correction 
Av/v 0 
(percent) 

Plate in direction of 

— > 0 

11.4 according to 

flow with finite 


Ferri 

thickness 


11.4 according to 
Lamia 

Elliptical cylinder 

0.1 

11.4 according to 
Ferri 

4.2 according to 
Lamia 


The numerical values found by Lamia, on account of the great 
simplification used, can only serve as guides for the approximate 
values of the corrections. The tunnel boundary conditions enter his 
computation only partly because equal flow was assumed through the 
cross section at the model position and the cross sections far ahead 
of and behind the model. This inexact taking into account of the 
effect of the tunnel wall led to the result that, in accordance with 
his computations for all Mach numbers, as for example for M = 0, 
the closed and the open tunnel have equal corrections of the velocity 
but with opposite sign. For compressible flow it is known, however, 
that the open tunnel requires only about one-half or one-fourth as 
large corrections as for the closed tunnel (reference 3, pp. 54-58). 
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2 . Approximate Computation of Frahke-Weinig 

A further approximation for the two-dimensional problem is given 
by Franke and Weinig (reference 4) . In their investigation they 
strictly consider the boundary conditions along the tunnel wall and 
take into account the compressibility in a manner similar to Prandtl. 
The velocity correction due to the displacement of the model is given 
in the following form: 


Av x _ f it 

v 0 6 h 2 (l~M 2 ) 3 / 2 

_ 1 
“ 3 

where 

F cross-sectional area of profile 

h distance of tunnel wall 


2^ V above + v below^ 


0 

m 


v upper v lower velocity at upper and lower tunnel wall, 

respectively, in plane of model 


v 0 


velocity far ahead of model 


M = v/a 


Mach number. 


tunnel velocity 
velocity of sound 


The second form of the velocity correction given previously, 
which refers to the disturbance velocity measured at the wall, is in 
many cases found to be very useful. This form is referred to also in 
the following investigation (see section IV) in cases in which the 
assumptions of the Prandtl principle in the neighborhood of the model 
are no longer satisfied although, with the aid of the wall velocity, 
a correction according to the Prandtl principle is still possible. 


3. Purpose of the Present Investigation 

The purpose of the present investigation is to give a useful 
approximation of the velocity and angle-of-attack correction, in 
particular for the circular- shaped closed wind tunnel. Because at 
high Mach numbers the wake produced by the model resistance greatly 
increases, an approximation to take account of this wake is developed. 
Further, the results obtained are checked with the aid of wind-tunnel 
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measurements so far as the available measurements permit such comparison 

In contrast to the work so far done, the computation is first 
carried out completely for the incompressible flow. The results 
thereby obtained are applied to the flow in compressible media by 
correlating each flow picture of the compressible flow with a definite 
flow picture in the incompressible flow with suitably modified tunnel 
and model dimensions so that for both flo.w pictures the corrections 
are either the same or stand in definite relation to each other. This 
computation method has the advantage of great clarity of computation, 
for the investigation is split up into two independent partial investi- 
gations. Furthermore, a large number of already existing correction 
computations for incompressible flows can, according to the same 
principles, be applied in a simple manner to compressible flows so that 
a good portion of the existing data is not lost but can be further 
utilized. 

The investigation was carried out in general form both for bodies 
of rotation and for wings of various spans. The numerical data given 
refer, however, predominantly to the case of the flow with rotational 
symmetry. In a further report that is soon to follow the numerical 
data will be supplemented. 


II. CORRECTION OF THE FLOW VELOCITY DUE TO MODEL 
IN CASE OF INCOMPRESSIBLE FLOW 

1. Equivalent Dipole or Source Strength of Model 

t 

(a) Equivalent dipole strength of model without wake. 

The flow about profiles or about bodies of arbitrary shape can, 
as is known, be simulated by definite arrangements of sources and 
sinks or of dipoles ( doublets ), the strength and distribution of which 
are so chosen that the resulting streamline coincides with the outline 
of the body investigated. For large distances from the dipole or 
source- sink system, the disturbance velocities produced by them can be 
shown to be equal to that of a single equivalent dipole that is located 
at the center of gravity of the system of singularities and whose 
strength is determined as follows- 1 -: 

Intensity of equivalent dipole 

M + = 2 2 (a Q) = 2 M 

^The equation given holds not only for doubly symmetrical but also for 
arbitrary body shapes . 
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where 

Q strength of source or sink 

2a distance between source and corresponding sink 
M moment of elementary dipole 

On the assumption that the dimensions of the model are small as 
compared with the diameter of the wind tunnel, the disturbance velocity 
of a model at the tunnel wall can be represented by the effect of the 
previously defined equivalent dipole M*. For large models, whose 
chord is comparable with the tunnel diameter, the simple estimate 
previously given of the equivalent dipole strength is insufficient; for 
this case a computation will later be given that permits estimating the 
deviation from the simple equivalent dipole. 

The equivalent dipole for various shapes of bodies was computed 
by Glauert as a function of the maximum thickness and the ratio of 
thickness to length of the body (reference 3) . In the present report, 
Glauert 's results are so modified that the equivalent dipole strength 
becomes a function primarily of the volume of the displacing body; the 
effect of the thickness ratio and the difference between two-dimensional 
and three-dimensional flow then remain very small for slender bodies. 

In this new representation the dipole intensity is given by 2 : 


2 For cylinders, Glauert gives the equivalent dipole strength as 

M + = \ b n/a d max 2 v 0 

(reference 3, p. 53) . The given value referred to the volume Ay 
is therefore connected with the A value of Glauert by the relation 

* 2 

x « a max b \ 

Ay = - A 

For the circular cylinder, A = 1 and therefore = 2. For bodies 
with rotational symmetry, 

M + = ^ \ dmax 3 v 0 

(reference 3, p. 59) and therefore 

A V ~ 4 ~Y~ K 

For the sphere, Glauert gives A = 1 so that Ay = 1.5. 



NACA TM 1300 


7 


= 2 2 (a Q) = Ay V Vq (1) 

where 

V volume of hody 
vq approach velocity 

For bodies with rotational symmetry the equivalent dipole is to be 
applied at the center of gravity of the body. For cylinders (t/b — * 0) 
the previously computed equivalent dipole strength is to be distributed 
along the axis of gravity of the cylinder so that along this axis there 
is a uniform dipole distribution dM* = db M + /b 

The factor Ay depends on the shape of the body and the thickness 
ratio. For very slender bodies' (d/t— »0) in two-dimensional and 
rotational symmetry cases, Ay has the same value Ay = 1. 

For bodies with elliptical outline, the factor Ay is given in 
figure 2 as a function of the thickness ratio for the two limiting 
cases of rotational symmetrical bodies and elliptical cylinder (t/b— *0). 
It is seen from the figure that, for slender bodies such as occur in 
airplane structures. Ay differs only by a slight amount from unity, 
for example. 


wing with d/t 

= 0.15 

Ay = 1.15 

fuselage with 

d/t = 0.25 

Ay = 1.065 


For the very thin plate transverse to the air stream, the volume 
of the body is zero. Because, however, such a plate will nevertheless 
give a displacement of the streamlines, the factor Ay must tend to 
infinity. If the computations in the literature 3 are used for the flow 
about a plate, the product is obtained for the plate with constant 
height h and very large span b (h/b — *0) : 



and for the circular disk with the diameter d: 



3 


For very wide plate: 


reference 5 . For circular plate : 


reference 6. 
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(b) Equivalent source strength to take into account the 
wake displacement behind a drag body 

Behind bodies which a flow drag connects with an energy loss 
(but not an induced drag) , as is known, a wake or dead-water region 
is formed that evidently increases the displacement of the body. 

(See fig. 3.) This dead-water region has its source in the neighbor- 
hood of the body and extends downstream to infinity, with increasing 
mixing with the normal flow. It is useful to represent the effect of 
this dead water for points at a great lateral distance away by a 
system of sources that arise at the place of the drag body and displace 
the normal stream behind the body by the same amount as the wake. 
Between the strength Q of this source system and the drag of the body, 
under certain simplifying assumptions, a relation can be given which 
may serve as an estimate for the displacement due to the wake^. 


^According to the momentum theorem, if the undisturbed pressure is 
assumed to prevail at the cross section considered behind the body, the 
drag W of a body is given by 



Av = mass flow per second X velocity loss 




= ~ P 



df 



(See fig. 3.) The integral contained in the last equation gives, 
however, precisely the additional strength per unit volume Q of the 
equivalent source that, on the assumption of potential flow, displaces 
the streamlines at a great lateral distance from the wake by the same 
amount as the wake, that is, 


Q = 




df or W = ~ p Vq Q 


By transformation 



w V Q 
P vq 2 


^ws v o / 2 


is obtained. The same relation between the source strength and the 
drag has already been given by reference 7 (p. 32) . 
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On the assumption that the static pressure in the wake differs 
only slightly from the pressure of the undisturbed flow, there is 
obtained, by neglecting small quadratic. terms, equivalent source 
strength: 

Q = \ v 0 f ws ( 2 ) 

where f wg = c w F _ represents the harmful drag of the body (without 
account taken of the induced drag) . This relation between the source 
strength and the drag area is independent of the shape of the body and 
holds both for two-dimensional and for three-dimensional flow. 

The exact location of this equivalent source is not uniquely 
determined. It appears admissible, however, to assume the source as 
located at the center of gravity of the volume . Although the position 
of the equivalent source may not be correctly given by this assumption, 
the error thus introduced is small provided that the equivalent is used 
only to represent the flow relations at a great distance from the source 
This assumption is, however, satisfied for obtaining the tunnel correct- 
ions because the model dimensions are always small in comparison with 
the tunnel diameter . 

In the derivation of the preceding equivalent source strength, the 
assumption was made that the static pressure in the wake differed only 
slightly from the pressure of the undisturbed stream. At a large 
distance behind the drag body this assumption holds quite well, so that 
to a very good approximation the previously given equivalent source 
strength remains the same although the dead-water region constantly 
expands. Immediately behind the wing, however, the static pressure 
may differ greatly from the undisturbed pressure of the flow so that 
deviations from the simple relation given above between the drag and 
the source strength are obtained. It. was shown by Muttray that the 
source strength decreases rapidly from a maximum value at the wing 
trailing edge to a final constant value behind the wing (reference 7, 
fig. 22) . This decrease means, however, that, to the previously 
computed equivalent source strength, a system of additional sources and 
sinks behind the wing is to be added. At a large distance from the 
wing, as for example at the tunnel wadi, this additional source system 
acts approximately like a single dipole. The action of this additional 
wake dipole is taken into account in the following correction compu- 
tations if the magnitude of the tunnel corrections from the disturbance 
velocities measured at the wall is determined according to equation (4). 
The error made in assuming that the location of this wake dipole does 
not coincide with the profile center of gravity is only of small 
significance on account of the effect of the wake corrections . 
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The same consideration applies also to the additional displacement 
due to the compression shocks, so that their effect can enter approxi- 
mately into the corrections provided the shock length is small as 
compared with the tunnel diameter. 


2. Disturbance Velocity at the Tunnel Center and at the Wall 
Due to the Model and the Dead Water 

Because, according to the assumption made, the model dimensions 
are small in comparison with the tunnel diameter, the disturbance 
velocities at the tunnel wall can be represented by the previously 
determined equivalent dipole and the equivalent source. The boundary 
condition which is to be satisfied by the presence of the tunnel wall 
for the closed tunnel is that all velocity components normal to the 
tunnel wall must vanish. 

This requirement is satisfied, as is known, by superimposing an 
auxiliary flow with a velocity field on the flow about the model in 
the free airstream so that the composite flow satisfies the previously 
given condition. The computaton for the circular tunnel is quite 
complicated. In the present case the solution was determined by a 
method given in reference 8 (p. 250) of developing the field of the 
radial velocities along the tunnel wall and the velocity field of the 
auxiliary flow into a Fourier series . By comparing the coefficients 
of the two Fourier series, the constants for the velocity field of 
the auxiliary flow were then obtained so that the velocity at 
arbitrary points of the flow field could be determined. As velocity 
correction, that velocity is taken which induces the auxiliary flow at 
the location of the model . 

The results of this computation are given in the present report 
for the additional velocities at the tunnel center and for the 
symmetrical case (with the equivalent dipole and equivalent source 
at the tunnel center) and for the additional velocity at the tunnel 
wall, because for these arrangements the computation is still rela- 
tively simple, only the first term of the Fourier series being taken 
into account. For several nonsymmetrical arrangements with respect 
to the excess velocity at the tunnel wall, numerical values are also 
given for which the first four Fourier terms entering the computation 
were taken into account. It is still necessary to check whether the 
further Fourier terms have an effect on the results^ . 


^A further report will be issued on . the required corrections of the 
given numerical data when the higher Fournier terms are taken into 
account in which the computation procedure will also be described in 
detail . 
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(a) Disturbance velocity due to the model. 

In the notation of Glauert, the velocity correction to be applied 
at the tunnel center is -written as follows: 



( 3 ) 


where 

V volume of body 

D tunnel diameter 

As has been previously explained, Xy represents a factor that 
takes into account the model shape. (See fig. 2.) The factor Ty 
depends on the tunnel shape and on the ratio of the model dimensions 
to the tunnel diameter. For a closed circular tunnel the factor Ty 
is given in table 1 for several typical cases. 


In the method given by Weinig, the additional velocity at the 
center of the tunnel may also be represented as a function of the 
additional velocity at the wall of the tunnel. As in the two- 
dimensional case computed by Weinig, it may also be shown for the 
closed circular tunnel that the disturbance velocities occurring at 
the plane of the model at the tunnel wall Ay^ due to the model 
bear a definite ratio to the disturbance velocity Av x at the tunnel 
center. (See fig. 4.) If, for the comparison velocity Av.^, the 
arithmetical mean of the velocities above and below the model is 
chosen, this mean value is independent of the lift because a potential 
vortex at the model at these two points induces equal but opposite 
circulation velocities . The wall velocity, according to figure 4, 
is composed of two parts : namely, a part that directly represents 

the disturbance velocity of the equivalent dipole and a part that is 
due to the wall effect. The required velocity Av x at the tunnel 
center may therefore be written as 


Av x 

v 0 


Av- 


m 


xw 


( 4 ) 


The factor m in this equation depends on the tunnel shape and on the 
ratio of the model dimensions to the tunnel diameter. For the closed 
circular tunnel the factor m for several typical cases is given in 
table 2 . 


12 


NACA TM 1300 


For comparison, it may be mentioned here that for a wing with infinite 
span between two walls, Weinig gives the value m = l/3. The factor m 
is likewise found to depend to a relatively large degree on the tunnel 
shape . 


(b) Correction for models of very large chord 

The numerical values given in the preceding section (a) for the 
velocity corrections were computed on the assumption that all the model 
dimensions except the wing span were small compared to the tunnel 
diameter. It is often necessary, however, to obtain an idea of how 
the data are modified if the preceding assumption no longer holds true 
with respect to the model chord. 

In order to investigate this effect on bodies of rotation in the 
center of the tunnel, the additional velocities at the tunnel center 
and at the wall were computed for various distances of the source-sink 
system replacing the body by the method given in section ( a) . The 
distribution of the sources and sinks for these bodies was chosen as 
shown in figure 5. With this arrangement, the ratio F p ro file/^ ^ 
was found to be equal to 0.75 as corresponds to the usual shapes in 
airplane structures and to wings and fuselages . Because for these 
bodies the magnitude of the equivalent dipole is known from equation (1) 
M 4 " = £ (2aQ), the disturbance velocities at the center of the tunnel 
and at the wall can likewise be immediately given by equations (3) and 
(4) if the effect of the large model chord is neglected. Through a 
comparison of the results, the numerical values for Ty and for the 
ratio m (equations (3) and (4)) were found to require a correction. 

The correction is represented in figure 6. For very large model chord, 
the additional velocities at the tunnel center are seen to differ only 
by a small amount from the value for vanishingly small models. The 
wall velcoity Av^ and therefore also the ratio m of wall velocity 
to center velocity depend, however, to a great extent on the model 
chord. For example, for a wing whose chord is equal to the tunnel 
diameter, the center velocity is 89 percent of the velocity obtained 
on the assumption of vanishingly small model chord, whereas the wall 
velocity is 66 percent of the corresponding value. 

The previous investigations for models of large chord refer only 
to bodies with rotational symmetry and to wings with vanishingly small 
span in comparison with the tunnel diameter. Strictly speaking, it 
would be necessary to carry out also corresponding computations for 
wings with various ratios of the span to the tunnel diameter. Before 
the results are extended in this respect, however, it may be assumed 
as a useful approximation that the correction values given in figure 6 
hold also for wings with finite ratio of span to tunnel diameter. 
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(c) Disturbance velocity due to the wake 

The effect on the tunnel wall of the_ dead-water region formed 
behind a drag body and therefore the effect on the tunnel correction 
can be given approximately by an equivalent source at the location of 
the model. (See section II, 1, (b).) If, for this flow, the boundary 
conditions at the tunnel wall are satisfied by superposition of the 
source flow in the free air stream on an additional velocity field 
according to section II, 2, (a), it is readily seen that in the closed 
tunnel at the model location no velocity component arises in the 
approaching flow direction because of the additional flow field; only 
for unsymmetrical arrangement of the source in the tunnel will a 
velocity component arise at right angles to the approaching flow 
direction, which leads to a change of angle of attack. This change of 
angle of attack is also, however, in general without significance, for 
models are practically never mounted unsymmetrically in the wind tunnel. 


Although no velocity component in the the flow direction is induced 
by the additional velocity field at the model location, a correction for 
the approach velocity in the closed tunnel is necessary. The relations 
are seen most simply with the aid of figure 7 in the case of the 
infinitely wide tunnel. The velocity field arising from the source 
in the tunnel and from the external sources reflected in the wall is 
characterized by the fact that the velocity components in the model 
plane in the flow direction vanish. At an infinite upstream and 
downstream distance, however, a parallel flow is formed with the 

velocity v=*I 0— For the flow in the tunnel at a large distance from 
2 F K 

the source this means that the approach velocity is no longer v' 0 but 

V/-, = v'n - — —■ The velocity vn therefore increases up to the plane of 
0 0 2 F k l q 

the model by the amount Av x = «■ — and up to a section very far behind 

r K 

the body by Av = Q/Fj^. The required approach- velocity correction in the 
closed tunnel is, therefore, with the air of equation (2), obtained as 


Av x = 1 Q = 1 £ws 

v 0 2 F K v 0 4 F IC 

where 

f ws = c w F harmful drag area 
Fj£ tunnel cross section 


( 5 ) 
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The distribution of this correction velocity over the model cross 
section is uniform; that is, at the tunnel center and at the wall there 
is the same velocity increment Av x . In the equation 

Av x/ v 0 = m Av xw/ v 0 
the factor m therefore has the value m = 1. 

This equation for the velocity correction due to the wake and 
magnitude of the factor m is independent of the shape of the closed 
tunnel and therefore holds for square as well as for circular tunnels, 
etc . 


In a. corresponding manner it may be shown simply that in wind 
tunnels with open test section no velocity correction is required to 
take into account the displacement due to the wake. The boundary 
condition of the free jet would be satisfied by reflecting sources 
and sinks alternately at the tunnel wall, in a similar manner to that 
indicated in figure 7. At the center of the tunnel no direct velocity 
component is obtained in the flow direction. At an infinite distance 
upstream and downstream of the plane of the source the velocity 
components due to 'the sources and the sinks reduce to zero as the 
velocities induced by them mutually cancel . 


III. VELOCITY CORRECTION DUE TO THE MODEL AND THE WAKE AT 
HIGH SUBSONIC VELOCITIES WITHIN THE REGION OF VALIDITY 
OF THE PRANDTL PRINCIPLE 

The equations given above for the tunnel correction hold only for 
the incompressible flow. With the aid of the Prandtl principle, 
however, a given tunnel with the model in a compressible flow can be 
correlated with another tunnel with models of modified dimensions in 
an incompressible flow. Between the velocities in the two correspond- 
ing tunnels definite relations hold so that, with the aid of these 
relations, the flow in the compressible medium can be reduced to the 
flow in the incompressible medium. The form of the Prandtl principle 
which is here used is as follows (reference 9) : 

At each point of the compressible flow the same potential holds 
as compared with the corresponding point of the incompressible flow 
and the same velocity in the flow direction (x-direction) with the 
v eloci ties at right angles to this direction reduced by the factor 
aJ±-V?- . The corresponding points in the two flow fields are eonnnected 
by the following relations : 
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x incompr ~ x compr 
yincompr = ycompr 


z lncompr — //l-M 2 

2 compr 


where 

M Mach number. 


tunnel velocity 
velocity of sound 


1 . Model Displacement 

For a wing with large aspect ratio B/t in a closed wind tunnel, 
the potential field for incompressible flow is shown in figure 8. The 
model is assumed to be slender enough and to possess a sufficiently 
sharp leading edge that the Prandtl principle is satisfied along the 
entire profile contour.. The profile contour is given by the boundary 
streamline of a definite source-sink arrangement so that the potential 
lines are continued to the interior of the surface enclosed by the 
boundary streamline . On the potential lines 0 = const so that they 
must intersect the tunnel wall and the profile surface at right angles, 
for at these bounding surfaces no normal velocity components exist . 


If this potential field is now transformed according to the 
Prandtl principle, that is, if the lines <f> - const are separ ated 
farther apart in the y and z directions by the factor l/ <^/l-M 2 , 
and the x component is kept constant, a new potential field is 
obtained for the compressible flow. (See fig. 8.) In this distorted 
potential field, condition for the tunnel wall, namely, that all 
potential lines must be at right angles to the flow direction, can be 
seen to be satisif ed along a new cont our whose y and z coordinates 
are increased by the factor 1 /a/ 1-M z . For circular tunnels this 
satisfaction means, for example, that the diameter of the tunnel in 
the compressible flow is greater than the comparison tunnel in the 
incompressible flow in the ratio 1/ l-M 2 . 


In the transformation according to the Prandtl principle the wing 
in the tunnel receives definite modifications. If, for the slender body, 
the assumption is made that, within the bounding potential line which 
replaces the profile, the potential lines may be replaced by their 
tangents at their points of intersection with the x-axis, it is seen 
tha t the inclination of this tangent becomes steeper in the ratio 
l/Vi-M 2 . This assumption means, however, that the stream lines in 
the neighborhood of the x-axis, as also in particular the bounding 
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stream, have flatter slopes dy/dx by the factor Vl-M 2 than before 
the distortion. The profile shape determined by the bounding stream- 
line is therefore thinner after the distortion in the ratio /y/l-M^ . 

The span of the wing with large aspect ratio B/t is determined by 
the length of the segment with sources and sinks or dipoles distributed 
on it. Th is di stance is likewise increased by the distortion in the 
ratio l/ iyl-M 2 so that the compari son w ing in the compressible flow 
has a greater span in the ratio l/ Vl-M 2 . 

Because the profile se ction has, through the transformation, become 
more slend er in the ratio Vi -M 2 but the span has become larger in the 
ratio l/ Vl-M 2 , the initial wing and the corresponding wing of the 
compressible flow have the same volume. It is also known for slender 
bodies of rotation that in the transformation according to the Prandtl 
principle their contour is, to a first approximation, unchanged so 
that in this case also the volume undergoes no change (reference 8) . 

Because in the Prandtl transformation the velocity components in 
the flow direction remain the same, a condition that holds also without 
restriction for the velocity components due to the model and the tunnel 
wall, the same corrections on the approach velocities are to be used 
for the bodies in the compressible flow as for a body transformed in 
the manner explained ab ove i n an incompressible flow in a tunnel 
reduced by the factor Vl-M 2 . 

If, therefore, a model that at small velocities has the approach- 
velocity correction Av x /vo for small velocities is placed in a closed 
tunnel and the Mach number continuously increased, the velocity 
correction according to equation (3) increases in the ratio 
l/(l-M 2 ) 3 / 2 , that is. 



(1-M 2 ) 3 / 2 


' Av x \ 

v 0 / incompr 


( 6 ) 


In the preceding relation, the factors Ty and Ay contained 
in equation (3) were assumed not to vary in the transformation accord- 
ing to the Prandtl principle. For Ay a small displacement actually 
occurs that is due to the change in the effective thickness ratio. 
According to figure 2, however. Ay depends only to a slight extent 
on the thickness ratio, so that the change in Ay may in general be 
neglected without great error. 


The factor Ty does not vary as long as the model chord is small 
in comparison with the tunnel diameter. With increasing Mach number, 
however, the effective ratio of model chord to tunnel diameter rises as 

l/ Vl-M 2 , because the model chord remain s unch anged while the tunnel 
diameter becomes smaller by the factor V 1-M 2 . When the Mach number is 
increased, the question is thus raised of the correction of models with 
large chord in a closed tunnel. According to figure 6, the correction 
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of the factor Ty can he estimated by letting the effective ratio of 
the model chord to the tunnel diameter increase by the factor l/ Vl-M 2 . 

As long as the model chord can still be considered small in compari 
son with the tunnel diameter, the factor m in equation (3) undergoes 
no change. In the transformation according to the Prandtl principle, 
the same ratio of model span to tunnel diameter occurs for the compari- 
son tunnel as for the initial tunnel, so that for both systems the same 
factor according to table 2 is to be used. If in the transformation 
too high ratios of chord to tunnel diameter result, however, then the 
value of m is to be corrected with the aid of figure 6. 


2 . Additional Velocity Due to the Wake 

The same consideration for the computation of the velocity correct- 
ion due to the model displacement can also be applied for the additional 
displacement due to the wake . In a purely formal manner, on the 
assumption that the drag surface does not vary with increasing tunnel 
velocity, there is obtained from equation (3) : 


(%) 


compr 


1-M C 


Ps) 


incompr 


( 7 ) 


It is thus seen that, in contrast to the cases so far given, the 
tunnel correction does not increase by l/(l-M 2 ) 3 / 2 but by l/(l-M 2 ). 
Because, in this computation, the effect of the wake was replaced by the 
flow about a source which, as is known, gives in a parallel stream a 
bounding stream line in the form of a "half-body," the preceding result 
means that for a body of small chord, for example, a fuselage in a 
tunnel, the velocity corrections increase by the factor l/(l-M 2 ) 3 / 2 j 
for a half-body in a closed tunnel the velocity corrections increase 
only by the factor l/(l-M 2 ). 




ss/////////////// 


Fuselage in 
tunnel 
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Ha If -body 
in tunnel 


The reason for this different behavior of the two body shapes is 
to be sought in the fact that for flow with rotational symmetry, the 
radial velocities acting on the wall at points a small distance upstream 
or downstream of the model for short bodies (equivalent dipole) decrease 
approximately in inverse ratio to the fourth power of the distance from 
the body, whereas for the half -body (equivalent source) they decrease 
only in the inverse ratio of the second power of the distance from the 
body. If, therefore, corresponding to the Prandtl principle, the tunnel 
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diameter in an incompressible flow is reduced for both bodies to the 
same extent, it is seen immediately that in the case of the h alf -body 

the wall effect must increase by a less amount than for the short body® 

/ 

If the velocity correction is again referred to the wall velocity, 
the factor m in equation (4) is still to be set equal to 1 for the 
region of validity of the Prandtl principle, because when the wake is 
taken into consideration, this factor is independent of the magnitude 
of the tunnel diameter. (See section II, 2, (c).) 


IV. VELOCITY CORRECTION DUE TO TEE MODEL AND THE WAKE FOR 
THE CASE WHERE THE PRANDTL PRINCIPLE ALONG THE MODEL 
SURFACE DOES NOT HOLD TRUE 

The assumption for the validity of the Prandtl principle is that 
the additional velocities due to the body are small compared with the 
tunnel velocity. In many cases, particularly at high Mach numbers, 
this assumption is not satisfied along the surface of the body. At 


%he radial velocities of source and dipole decrease according to these 
considerations, for small values of x, by the second or fourth power 
of R; the addition al vel ocities increase, however, with the second or 
third power of l/ ^ 1-M 2 . This difference in the exponential relation 
is to be ascribed to the fact that at a great distance ahead of or 
behind the source and the dipole, the radial velocities in each case 
vary by the same amount so that an equalization of the power exponent 
occurs . A simple check computation for the amount of the velocity 
increase is possible for the half -body. If a half -body is in a closed 
tunnel, the stream is parallel far ahead of and behind its nose. The 
velocity increase is computed simply by the decrease in the cross 
section df. From the known equations 

dp _ k M 2 df 


and 


P 




in agreement with the relation previously found from the Prandtl 
principle : 


dv 


1 df 
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some distance ahead of or behind the body the additional velocities 
have become small enough, however, that for the flow outside of a 
definite boundary line the Prandtl principle may again be applied. 

(See references 10 and 4.) In figure 9, for example, the region , 
enclosed by the boundary line no longer follows the Prandtl principles 
so that the transformation law is unknown for the streamlines within 
the boundary and for the body contour. 

This fact is now applied to the flow in a tunnel, first for the 
simplest case of the flow in an infinitely wide tunnel with parallel 
walls . In an incompressible flow the boundary conditions at the 
tunnel wall are satisfied for the dipole replacing the body by an 
infinite reflection of the equivalent dipole at the tunnel wall . In 
the neighborhood of the dipole the disturbance velocities are so large, 
however, that within a certain limiting region the Prandtl principle 
is no longer applicable (fig. 10) . The transformation of this 
incompressible flow field according to the Prandtl principle is possible, 
however, for the regions lying outside with boundary lines; within 
the boundary lines, as stated above, a transformation must be made 
according to an as yet unknown law so that the body contour for the 
compressible flow is likewise unknown. According to the Prandtl 
principle, the strength of the equivalent dipole can no longer be 
determined by computation for a given body. This gap can however 
be closed for the wind-tunnel test by a simple measurement. Because 
the additional velocity at the tunnel wall above and below the model 
is a measure of the effective dipole strength and the tunnel walls, 
moreover, lie within the region of validity of the Prandtl principle, 
the effective dipole strength can be determined by measuring the 
velocity at the wall. The required tunnel- correction velocity at the 
tunnel center, is as is known, the sum of the induced velocities of 
all reflected dipoles whose intensity is k n own from the measurement 
of the wall velocity and for which the tunnel center belongs to the 
region of validity of the Prandtl principle. This knowledge means, 
however, that between the wall velocity and the approach- velocity 
correction at the tunnel center the same relation holds as in an 
unrestricted region of validity of the Prandtl principle, that is, 


Av x/v 0 = m Av^/vq 

This relation is not restricted to only the flow between two 
walls. As can be readily seen, it can be applied also for the closed 
tunnel of arbitrary cross-sectional shape, for example, for the 
circular tunnel. The only difference consists in the fact that in the 
flow between two walls the additional potential to satisfy the tunnel 
wall conditions is produced by the reflected singularities, whereas 
in the most general case the additional potential arises also from the. 
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singularities outside of the tunnel; the strength and the location of 
these singularities cannot be found, however, by simple reflection. 

As long as the model dimensions are small as compared with the 
tunnel diameter the values to be assigned to m are those given in 
table 2. For models with large chords and for high Mach numbers, 
however, an estimate for the correction of the factors m due to very 
large ratio of the chord of the model to the tunnel diameter must be 
made with the aid of figure 6. 

The relations given between the velocity at the wall and the 
tunnel- correction velocity for flows about models, for which in the 
neighborhood of the model the Prandtl principle no longer holds, can 
also be derived for the velocities induced by the wake. 

The previously described extension of the tunnel-correction 
computations to Mach numbers for which the Prandtl principle in the 
neighborhood of the model no longer holds is naturally inapplicable 
without restriction up to a Mach number M = 1. An upper limit of the 
applicability is given, for example, if the sound velocity is attained 
or exceeded at the tunnel wall . In this limiting case the assumptions 
of the Prandtl principle are no longer satisfied in the neighborhood 
of the tunnel wall so that the entire transformation process is 
inapplicable- How closely this upper limit may be approached in the 
wind tunnel measurement without fundamentally altering the pressure 
field about the profile can be determined by experiment. This limiting 
Mach number will depend, among other factors, on the degree of 
obstruction of the tunnel by the model and on the angle of attack of 
the model . Corresponding tests to determine admissible model dimensions 
for definite limiting Mach numbers and angles of attack are at present 
being conducted at the DVL. 


V. TUNNEL CORRECTIONS DUE TO LIFT AT HIGH SUBSONIC VELOCITIES 
1 . Wing in Free Air Stream 

In order to consider briefly the flow relations about a wing with 
lift in a compressible medium, the wing will first be considered replaced 
by a vortex filament in a free air stream. The flow field about a 
potential vortex in an incompressible flow is schematically represented 
in figure 11 . If the velocities along the closed bounding line ABCD 
are added, the following integral represents the circulation of the 
vortex : 

n 


v ds = Z = circulation 
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The magnitude of the circulation is, as is known, independent of 
the path of integration as long as the vortex lies within the boundary 
line. In a medium of constant density there is then obtained for the 
lift per unit length 


A = p Z v 0 

where p is the air density. 

If the flow field of this vortex is now transformed according to 
the Rrandtl principle where there is again excluded a region in. the 
immediate neighborhood of the vortex center, the control surface ABCD 
is extended in the direction of the y axis. (See fig. 11.) At the 
same time, however, the veloc ity co mponent v y is reduced at corres- 
ponding points in the ratio 1-M^” so that the product of the path 
element by the velocity at the normal boundary lines undergoes no 
change as a result of the transformation. Because, however, in the 
x direction neither the lengths nor the velocities have changed, the 
circulation integral along the lines ABCD or A'B'C'D' remains 
unchanged. If the line A'B'C’D' is taken at a very large distance 
from the vortex, where the disturbance velocities and therefore the 
pressure and density changes in the flow have been reduced to vanishingly 
small values, the compressible flow in the neighborhood of the control 
line completely resembles an Incompressible flow, which means that for 
the lift per unit length the familiar relation again holds : 

A = p v 0 Z 

where p is the density of the medium at a great distance from the 
vortex. 

From the preceding consideration, the conclusion can therefore 
be drawn that in the Prandtl transformation the circulation and the 
lift referred to the span element are unchanged. 

For compressible flow the same lift per unit spa n is, as known, 
attained at an angle of attack reduced in the ratio A^/l-M^ . This is 
seen from the fact that in the transformation the velocity components 
normal to the flow direction become smaller while the components in 
the flow direction remain unchanged. The slope of the stream lines 
and therefore also the angle of attack thus become smaller In the 
same ratio . 

The same considerations may be applied also for a wing of finite 
span with elliptical lift distribution. A wing in a compressible flow 
is then to be compared with a wing in an incompressib le fl ow for which 
the span of the second wing is reduced in the ratio Vl-M 2 . At each 
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section at corresponding distances from the center of the wing the 
circulation and the lift per unit span element then agree for both 
wings with elliptical, lift distribution, whereas in the incompressible 
flow the angle of attack, b oth t he geometrical and the induced, are 
increased in the ratio 1 / <|/l-M^ . The comparison wing in the 
incompressible flow then has an aspect ratio (span b/chord t) imparied 
by the factor \/l-M2 . For the induced drag in the compressible flow 
the equation 

c a 2 F 

c wi - c a Aai - — - 
nevertheless still holds where 


F wing area in compressible flow 

b span in compressible flow 


because for the comparison wing ip the inc ompre ssible flow the induced 
angle of attack Act^ is increased by l/ ^/l-M - ^ while through the 
transformati on to a compressible flow this angle of attack is again 
reduced by fy/l-M^ . The two effects thus mutually cancel . 


The equations for the induced angle of attack 


Aojl 


fa _F_ 
n b 2 


and the induced drag 


-wx 


71 ^ 


retain their validity also for compressible flow. 


2 . Wing with Lift in the Wind Tunnel 

As has already been stated in corresponding cases, a wing in the 
wind tunnel with compressible flow can be made to correspond to a 
definite wing with equal lift in a similar reduced wind tunnel with 
incompressible flow. Between the correction velocities of the compari- 
son tunnel and the initial one the relation then exists that the Vy 
c orrec tions in the compressible flow must be made smaller in the ratio 
jJ l-M^ . In this manner all correction computations can be carried over 
to the compressible flow. 
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The results of the existing wind-tunnel-correction computations 
are, in general, represented in the following form: 


Angle of attack correction: Aa = 8 -g- 

where 

5 correction factor, function of ratio wing span to tunnel diameter, 
tunnel shape, and lift distribution 

F wing area 

Fr tunnel cross section 

In the Prandtl transformation the value of the factor 8 is not 
changed because the ratio of span to tunnel diameter, which determines 
its value, remains the same and so do the tunnel shape and the lift 
distribution. The wing area, however, becomes smaller in the ratio 
/y/l-M 2 , the tunnel cross-sectional area in the ratio (l-M 2 ) so that 
for the comparison tunnel in the in compr essible flow the angle of 
attack is greater in the ratio l/ y.l-M 2 . This increase is, however, 
again canceled in converting to the compressible flow because all v y 
velocities and angl es of attack become smaller in the Prandtl trans- 
formation by Vl-M 2 . 

The important result is thus found that for elliptical lift 
distribution the angle- of- attack corrections and therefore also the 
corrections of the induced drag can be dealt with in the same way as 
for incompressible flow. 

The assumption underlying the preceding general result is, 
however, that the dimensions of the wing with the exception of the 
span are small compared with the tunnel diameter and only the corrections 
at the location of the wing are considered. If such is not the case, 
as, for example, in the correction of a wing with large chord due to 
the stream curvature in the wind tunnel or in the computation of the 
downwash behind a wing, it is again advisable to make use of the idea 
of a comparison tunnel with its corrections determined. It is thus 
found, for example, that the correction due t o the stream curvature 
for large chord wings increases as l/yi-M 2 , the profile chord remain- 
ing the same in the Prandtl transformation, while the tunnel diameter 
becomes smaller by the factor f\/ l-M 2 . With increasing Mach number the 
ratios of the wing chord to the tunnel diameter also increase in the 
comparison. The corrections due to the flow curvature may therefore 
at high Mach numbers be of significance, although for an incompressible 
flow they may be entirely negligible. 
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VI. APPLICATION OF TEE WIND TUNNEL CORRECTIONS DERIVED PREVIOUSLY 
1. Superpos ability of the Individual Corrections 

In the present work the different factors for the wind tunnel 
corrections, like model displacement, wake displacement, and lift, 
were treated as though only one of these magnitudes was alone effective 
(for example, a displacement body without lift and drag, or a lifting 
vortex without displacement, and so forth) . In the wind tunnel test, 
however, the various factors enter in general together so that the 
question arises of the superposability of the individual corrections . 

In the incompressible flow the question can immediately be 
answered in that the individual factors are simply to be superposed 
linearly. Each individual correction can, as is known, be computed 
from the corresponding potential field. Because for incompressible 
flow the potentials can be superposed linearly, the corrections derived 
from them can similarly be superposed linearly. 

In the case of the compressible flow, a flow picture in the 
incompressible flow was made to correspond to each flow picture, with 
the aid of the Prandtl principle . Between the velocities and angles 
of attack of the two flow fields the familiar relations of the Prandtl 
principle hold. In the incompressible comparison flow field, the 
individual factors and the corrections may again be linearly superposed. 
Because the required corrections in the compressible flow differ from 
these comparison corrections only by a factor, the law of linear 
superposition of the corrections holds also for the compressible flow. 


2 . Application of the Corrections in the Wind-Tunnel Test 

(a) The corrections due to the lift are obtained according to 
the known equations from the measured lift. The corrections do not 
increase, for equal lift coefficient, with increasing Mach number as 
long as the model chord in the comparison tunnel remains small as 
compared with the tunnel diameter. The corrections due to flow 
curvature for large-chord mod els i ncrease, however, with increasing 
Mach number in the ratio 1 / \/l-M 2 so that these corrections can 
become of significance although for the incompressible flow they are 
entirely negligible. 

(b) The corrections due to the measured drag are computed 
according to equations (5) and (7). They increase for equal drag area 
fws = c w F as l/( 1-M 2 ) . 
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(c) The corrections due to the model displacement are, for small 
Mach numbers, estimated according to equations (3) and (6), as long as 
it is assumed that along the body contour the assumptions of the Prandtl 
principle are satisfied with sufficient accuracy. The corrections 
increase as 1/(1 -m 2)3/2_ if at high Mach numbers in the neighborhood 
of the model the Prandtl principle no longer holds, the correction 
velocity at the tunnel center is obtained from the additional velocity . 
measured at the wall above and below the model with the aid of equa- 
tion (4) and the factor m given in table 2. From the additional 
velocity measured at the wall, the part due to the drag is to be 
subtracted as this correction is already taken into account. In 
addition, the effect of the mounting is naturally to be taken into 
account as is done most simply by a calibration measurement. 

At high Mach numbers the previously given corrections depend on 
measurements that are obtained at a large distance from the model ' and 
that there have the same values as in the neighborhood of the model . 
Because the assumptions for the Prandtl principle used in these correct- 
ions are well satisfied at a large distance from the model, the previous 
corrections are also admissible when the Prandtl principle in the 
neighborhood of the model no longer holds . 


3. Sample Computations and Comparison with the Approximate 
Computations of Ferri and Lamia 

For the velocity correction due to the model displacement the 
previously derived equations give higher approximations than those of 
Ferri or of Lamia. In order to obtain an idea of the admissibility 
of the assumptions of Ferri or Lamia it is therefore convenient to 
compare, with the aid of a few examples, the various approximate 
computations. (See table 3.) 

From the comparison of the values in the table it is seen that, 
in accordance with expectation, the wind-tunnel corrections of Ferri 
are greatly overestimated. The values of Lamia, particularly for the 
tunnel with open section, also still lie considerably higher than the 
more accurate values of Franke-Weinig or Glauert-Gothert . 

The comparison of the closed circular tunnel with the closed 
tunnel with plane walls shows that in spite of the same ratios of 
model thickness to tunnel height the wind-tunnel corrections consider- 
ably deviate from each other. 
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4. Relations from the Adiabatic Equation as an Aid to 
Wind-Tunnel Corrections 

For the correction of dynamic pressure and Mach number, several 
of the relations obtained from the adiabatic equations are of import- 
ance . These are : 

Correction of the dynamic pressure: 


dq/q = (2-M 2 ) dv/v = - l/2 (2-M 2 ) dp/q 

Correction of the Mach number: 

dM/M = (1+ M 2 ) dv/v = - l/2 (1+ ^ M 2 ) dp/q 


where 

p static pressure » 

q pv 2 /2, dynamic pressure 

v velocity 
k 1.405 for air 

k-1 2 . 

The value of the expression in parenthesis (1 + — — M ) is repre- 
sented in figure 12. The factors \/l-M 2 , 1-M 2 and (l-M 2 )^/ 2 , which 
occur frequently, are also plotted in figure 13 . 


VII. COMPARISON OF THE COMPUTED CORRECTIONS WITH 
WIND TUNNEL CORRECTIONS 

To check the tunnel corrections determined above for high subsonic 
velocities, the high-speed DVL wind tunnel was available. This wind 
tunnel has a closed measuring section of 2.7-meter diameter and attains 
at about 50 percent of the available driving power the velocity of 
sound in the measuring section. 


1. Changes in the Wall Pressure through the Mounting of 

Rectangular Wings of Various Chords 

In order to learn the effect of the tunnel wall on the measure- 
ment values at high subsonic velocities, four rectangular wings of 
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equal profile, KACA 0015-64, with various chords were investigated in 
this tunnel,- the chords were t = 350, 500, 700 and 1000 millimeters. 

The span (B = 1.35 m) and the mounting were the same for all four wings 
investigated. The additional velocities produced by the tunnel walls 
cannot be directly measured. At most the increase in the surface 
pressure at various Mach numbers could be compared insofar as the change 
in the surface pressures in the free air stream at high subsonic 
velocities could be considered as known from some computations. In 
this case too, however, the observed increase in the surface pressure 
is due partly to the increase of the flow about the profile as a result 
of the compressibility effect and partly to the effect of the tunnel 
wall, so that the part due to the wall effect, obtained through split- 
ting the measured values, is at least uncertain. 

A useful method of checking the corrections is offered, however, 
by the measurement of the wall pressure in the model plane. The 
computation of the wall pressure depends on the same assumptions as the 
computation of the correction velocities at the tunnel center, as has 
already been explained in the previous sections. This close connection 
■ expresses itself also in the fact that the pressure changes due to the 
model measured at the wall stand in qufte definite relation to the 
correction velocities at the center of the tunnel. (See equation 4.) 

In addition, the pressure changes at the tunnel wall are always greater 
than the pressure changes entering the correction computation so that 
the wall pressures are more easily susceptible to an accurate measure- 
ment. The tests were conducted by measuring the wall pressures pj. 
Pm, and p^- for various dynamic pressures. (See fig. 14.) The 
measuring stations for the pressures pj and Pm were uniformly 
distributed over the entire circular cross-section. For the measure- 
ment of p w in the test section, three close-lying holes were bored - 
in the wall in the model plane above and below the model and these were 
combined to give the arithmetic mean. It may be shown that the dynamic 
pressure and the Mach number in the plane of the model and the wall 
pressure p w without the model were functions only of the ratio 
(Pp - Plll)/Pi- These relations were determined by tests and plotted 
as calibration curves . The values determined from this ratio for the 
dynamic pressure and the Mach number are, in what follows, denoted as 
the uncorrected test- section values. If a model is mounted in the test 
section, the indicated wall pressure p w changes for equal pressure 
ratio (pp - Pixx)/pi- From this change of the wall pressure, the 
corrections to be applied can be determined. 

The wall-pressure changes measured for wings of various chords 
are plotted in figure 15 as functions of the corrected Mach number in. 
the test-section center. In order to eliminate the effect of the 
mounting, the wall-pressure change was not referred to the wall pressure 
of the free test section but to the wall pressure for a mounted wing 
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of 350-millimeter chord. In order to obtain a computational comparison 
the Prandtl principle was assumed to be valid also in the neighborhood 
of the model. Although this assumption does not hold with certainty 
if compression shocks occur at the model, that is, at a Mach n umb er 
above approximately 0.76, it may nevertheless be assumed that the 
admissibility of the assumptions can be checked from the trend of the 
computed and measured curves . This assumption holds at least for the 
greater part of the curve, which is below the critical Mach number of 
0.76. The computation of the wall-pressure changes was then carried 
out with the aid of equations (3), (5), (6), and (7) and tables 1 and 2, 
first without taking account of the measured wing drag and again with 
the drag taken into account. It was found that for the wings investi- 
gated, the wing drag becomes of significance only if it has increased 
greatly because of the compression shock. Nevertheless it even then 
has a small effect so that the error in taking account of the wing drag 
has not too great an effect on the curves . 

The effective Mach number with which the rise in the wall pressure 
is to be computed is not unique, particularly at high Mach numbers. 

The Prandtl principle requires that the mean Mach number of the flow 
field under consideration be substituted. In order to show the effect 
of this uncertainty, two curves were computed, one for the corrected 
Mach numbers at the tunnel center, which in any case had 'to be considered 
as too small in the neighborhood of the model. For the other curve the 
Mach numbers computed at the wall were used, which in the neighborhood 
of the body came close to the mean Mach number. From the curves in 
figure 15 it is seen that the computed curves on the whole agree well 
with the measurement points. The existing deviations all lie within 
a scatter range that can be explained by an inexact estimate of the 
effective Mach number . 

It may therefore be said in conclusion that the computations 
throughout agree with the measurements to the degree that may be 
expected from the assumptions made . 

With regard to the previously mentioned uncertainty in the deter- 
mination of the effective Mach number, it is further to be added that 
this uncertainty does not exist in the correction of 'the dynamic 
pressure and the Mach number. The principal factor of importance, 
namely, the ratio of measured wall-pressure change to the tunnel cor- 
rection, is independent of the corresponding Mach number. Only in the 
determination of the correction factors m corr /m and t v corr / v at 

large model chord and in the taking into account of the wake displace- 
ment does the Mach number enter . The deviation of the computation 
through uncertainties in these corrections should not be of great 
signif icance .' 
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2 . Effect of the Lift on the Wall Pressure 

The measurements described above on the rectangular wing were 
carried out for a symmetrical flow about the wing, that is, for zero 
lift . The fact that a change in the angle of attack to a first approxi- 
mation has no effect on the mean wall pressure is to be expected from 
the considerations of section II, 2. The circulation due to the lift 
produces on the upper and lower sides equal and opposite velocity and 
pressure changes so that the disturbance velocity through the circu- 
lation in forming the mean of the wall pressures on the top and under 
sides drops out. There is only an effect due to changes in the wing 
drag as a result of different angle of attack, which, however, because 
of the small effect of the drag on the wind-tunnel corrections, cannot 
be of great significance. 

These facts could be confirmed by tests. Figure 16 shows, for a 
rectangular wing of 500-millimeter chord, the wall pressure j? w as a 
function of the pressure difference pj - pm, already discussed 
previously for 0° and 5° angles of attack. It is seen that there are 
no systematic deviations between the measuring points of the various 
angles of attack that could not be ascribed to errors in measurement. The 
5° angle of attack means, however, that at high Mach numbers, in general, 
the limit of the angle of attack is of greatest interest in the wind 
tunnel experiment. This independence of the angle of attack means a 
considerable simplification in evaluating the results because the cali- 
bration curves for dynamic pressure and Mach number need not be corrected 
for each investigated angle of attack. 


3. Pressure Drop in the Test Section Due to Wake Displacement 

In order to check equations (5) and (7) for taking into account 
the wake behind the drag bodies, use was made of the fact, discussed 
in section II, 2, c, that at a large distance behind the model the 
disturbance velocity due to the wake is just twice as great as the 
corrections at the wing location according to the equation (5). It 
was further assumed that at the end of the test section, that is, about 
two wing chords behind the investigated wing of 500-millimeter chord, 
this final value is already practically attained. Under these assump- 
tions the additional pressure drop due to the wake as compared with the 
test section without obstacles could be computationally estimated. For 
the effective Mach number there were again substituted two values which 
correspond respectively to the Mach number in the' test section and the 
Mach number at the end of the test section. 

For this comparison too (fig. 17), the measuring values lie within 
the region that is described by the computed curves . The residual 
deviations can be ascribed to measuring errors or deviations in the 
effective Mach number. 
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VIII . SUMMARY 

(l) For wings with finite ratio of span to tunnel diameter an d 
for bodies of rotation in closed circular tunnels, the wind-tu nn el 
corrections due to the model displacement in an incompressible flow 
are computed. 

The corrections, in contrast to Glauert's method, are given as a 
function of the volume of the displacing body. In this method the 
effect of the contour shape for slender body shapes becomes vanishingly 
small,- in particular in the limiting case of very slender bodies the 
same form factor is obtained for three-dimensional and two-dimensional 
flows . 


(2) For incompressible flow the additional velocities due to the 
dead-water region behind the resistance bodies is represented by a 
simple equation. 


(3) With the aid of the Prandtl principle it is shown that for 
compressible flow the tunnel corrections due to the model displacement 
increase as l/(l-M 2 ) 3 / 2 and due to the wake for equal drag area as 
l/(l-M 2 ). The corrections due to the lift remain, for equal lift 
coefficient, unchanged provided the wing chord is small compared to 
the tunnel height. 


(4) On increasing the Mach number the corrections to take accou nt 
of the stream curvature for models with large chord rise as i/a/1-m 2 
so that these corrections have significance at high Mach numbers even 
though they are negligible at small Mach numbers . 


(5) The derived corrections for high Mach numbers remain valid 
if along the model surface the assumptions of the Prandtl principle 
are no longer satisfied. The limiting Mach number up to which the 
method is applicable is to be determined by wind-tunnel tests. 


(6) A comparison between computation and measurement shows good 
agreement as far as may be expected from the assumptions made. 


Translated by S. Reiss 
Rational Advisory Committee 
for Aeronautics 
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TABLE I - FACTOR T y FOR THE CLOSED CIRCULAR TUNNEL 1 (EQUATION 3) 
[b, span of rectangular wing; D, tunnel diameter] 


Model shape 

Body of 
rotation 

Rectangular wing 

B/D = 0 

B/D = 0.25 

B/D = 0.50 

B/D~= 0.75 

Factor T y 

1.02 

1.02 

1.04 

1.06 

1.10 


■''Glauert (reference 3, p. 58) gives for the factor T y after suitable 
conversion for bodies with rotational symmetry: 


T y = 0.797 X 4/rt - 1.016 
for closed circular tunnel, 

T y = -0.206 X 4/jt = -0.263 

for open circular tunnel. The deviations of these values from those 
in table 1 lie within the accuracy of computation. 


TABLE II - FACTOR m FOR THE CLOSED CIRCULAR TUNNEL (EQUATION 4) 


[B, span of rectangular wing; D, tunnel diameter] • 


Model shape 

Body of 
rotation 

Rectangular wing 

b/d = 0 

B/D = 0.25 

B/D = 0.50 

Factor m 

0.45 

0.45 

0.46 

0.49 
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TABLE III - VELOCITY CORRECTION IN VARIOUS WIND TUNNELS FOR 
A WING WITH ELLIPTICAL CROSS-SECTION WITH d/t = 0.10 
FOR A MACH NUMBER OF 0.75 


Tunnel shape 

Wing thickness 

Span 

Velocity correction 


Tunnel height 
(percent) 

Tunnel width 

(percent) 

Closed tunnel 
with parallel 
walls ; 

h = constant; 
b oc 

4 

1 

11.4 -- Ferri 

4.2 -- Lamia 

2.5 -- Franke- 

Weinig 
and Glauert- 
Gothert 

Open free jet 
with parallel 
walls; 

h = constant; 

b — » eo 

4 

1 

-4.2 -- Lamia 

-1.3 — Glauert- 

Gflthert 


Closed circu- 
lar tunnel 


4 


0.25 


1.2 


— Gather t 
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Figure 


Figure 2. 



.. - Wind tunnel' with model. One-dimensional velocity distribution according 

to Ferri . 
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(recomputed from reference 3). 
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Figure 3. - Model with vortex wake. 
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Figure 4. - Velocity distribution for a model in the wind tunnel. 



Figure 5. - Distribution of the sources and sinks for the hounding 
streamlines under consideration. 
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Figure 6. - Effect of the model chord on the factors Ty (equation (3)) 

and m (equation (4)) for source-sink body with ratio F/d X t = ~ 0.75 
(rotational symmetry or wing with ratio of span/tunnel diameter -> 0) . 
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■Dipole 
in the 
tunnel 



Figure 10. - Dipole in tunnel with "boundary lines for the case of the 

validity of the Prandtl principle. 




x 


. - Potential vortex in the parallel flow. 


Figure 11 





Mach number, M = v/a 


1-M 2 , and 


(1-M 2 ) 3 / 2 


as a function of the Mach number 
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Figure 14. - Wall pressures for determining the dynamic pressure 




Bill 


= pressure difference at the 
walls for wing of 350 mm 
chord for equal mass flow 
of air through test section. 
pv 2 /2 = dynamic pressure at 
center of tunnel (corrected) 


Test points: 

a Wing with 1000 mm chord ■ 
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Computed curves: 
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a = computed for Mach number 
at tunnel center 
b = computed for Mach number 
at tunnel wall 
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Te3t ■points: 


0 Without model 
• With model 




Calculated with 

corrected Mach 

numher a - j 

1 1 Tjk 

Test section with 
wing (t = 500 




Calculated with 
uncorrected 'Mac. 
numher 




LSr = je = : Test section without 
installations 

(test-mean curve) 

Mach numher (uncorrected) 
at center of test section , 


Figure 17. - 


Pressure drop Ap m the test section with and without the model 
Ap = Pj n it " Pfin = wall-pressure difference between the begin- 
ning and end of the test section. 

q-ad = p£ 0 £ - pij v j = adiabatic dynamic pressure at the center of 
the test section. 








